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Abstract 

This  paper  describes  calculations  of  the  electrical  losses  in  patterned  series-connected  solid  oxide  fuel  cells  (SOFCs),  i.e.  the  “segmented-in- 
series”  configuration.  Losses  due  to  cell  resistances,  electrode  ohmic  resistances,  interconnect  resistance,  and  shunting  by  a  weakly-conductive 
support  material  were  considered.  For  any  given  set  of  cell  dimensions  and  characteristics,  power  density  was  maximized  at  an  optimal  cell 
length  -  the  power  decreased  at  larger  cell  lengths  due  to  electrode  lateral  resistance  loss  and  decreased  at  smaller  cell  lengths  due  to  a 
decreasing  fraction  of  cell  active  area.  For  cell  lengths  well  above  the  optimal  value,  electrolyte  current  was  often  confined  to  the  portion  of 
the  cell  nearest  to  where  the  less  conductive  electrode  (i.e.  the  cathode)  connected  to  the  interconnect.  Assuming  dimensions  expected  for 
screen  printing,  i.e.  f»20  pm  thick  electrodes  and  lateral  print  accuracy  of  «T00  pm,  and  area  specific  resistance  values  typical  of  SOFCs, 
optimal  cell  lengths  typically  ranged  from  1  to  3  mm.  Shunting  currents  increased  with  decreasing  cell  lengths,  but  were  found  to  have  little 
effect  on  power  density  assuming  partially-stabilized  zirconia  supports  and  temperatures  <  800  °C. 
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1.  Introduction 

The  “segmented-in-series”  solid  oxide  fuel  cell  (SS- 
SOFC)  design  has  often  been  implemented  in  a  tubular  con¬ 
figuration,  with  the  cells  delineated  in  bands  around  the 
tube,  leading  to  the  term  “banded  tubular”  [1].  Fig.  1  shows 
schematically  a  cross-sectional  view  of  this  geometry.  Re¬ 
cently,  a  similar  design  using  flattened  tubes  instead  of  circu¬ 
lar  tubes  has  also  been  reported  [2].  Cell  lengths  were  typi¬ 
cally  >  10  mm  in  both  circular  and  flattened  tube  designs.  In 
addition  to  the  general  advantages  of  tubes  -  seals  are  not 
required,  and  the  interconnect  is  an  integral  layer  rather  than 
a  separate  piece  -  flattened  tubes  allow  for  higher  packing 
density,  thereby  yielding  higher  power-to-volume  ratios  and 
facilitating  the  use  of  low-cost  planar  deposition  methods 
such  as  screen  printing. 

General  predictions  regarding  the  performance  of  SS- 
SOFCs  have  been  made  based  on  modeling  of  the  electri¬ 
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cal,  electrochemical,  heat  and  mass  transfer  processes  [3].  In 
a  recent  paper,  we  suggested  that  SS-SOFCs  could  provide 
improved  performance  if  the  cell  lengths  were  decreased  be¬ 
low  10  mm,  and  demonstrated  the  fabrication  of  1-2  mm  long 
cells  by  screen  printing  [4].  A  simple  calculation  also  showed 
that  the  increased  number  of  interconnects,  due  to  the  larger 
number  of  cells  per  unit  length,  should  not  be  detrimental 
to  performance.  An  estimate  of  the  support  shunting  current 
was  also  given.  However,  neither  of  these  papers  considered 
in  detail  the  limitations  on  SS-SOFC  geometries  imposed 
by  practical  processing  methods,  or  made  predictions  con¬ 
cerning  the  optimal  geometries  and  their  expected  perfor¬ 
mance. 

In  this  paper,  we  present  quantitative  performance  predic¬ 
tions  for  SS-SOFCs  as  a  function  of  cell  and  interconnect 
geometry,  support  material,  cell  area-specific  resistance,  and 
electrode  sheet  resistance.  Most  of  the  calculations  were  done 
assuming  ceramic  layer  fabrication  by  screen  printing,  where 
the  minimum  feature  size  is  « 100  p.m.  The  results  provide 
guidance  regarding  optimal  SS-SOFC  geometries  and  their 
expected  power  densities. 
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2.  Calculation  methods 

Two  main  calculations  were  carried  out.  First,  the  elec¬ 
trode  resistance  losses  across  segmented-in-series  cells  and 
interconnects  were  calculated.  Second,  the  shunting  current 
between  cells  across  the  support  was  calculated.  Cell  dimen¬ 
sions  were  chosen  to  reflect  the  capabilities  of  screen  printing, 
a  commonly  used  SOFC  fabrication  method.  In  terms  of  lat¬ 
eral  print  resolution,  all  lines  and  spaces  were  chosen  to  be  at 
least  0.1  mm  to  reflect  what  is  achievable  with  typical  screen 
mesh  counts  (230  and  325),  as  well  as  visual  pattern  align¬ 
ment.  Layer  thicknesses  were  also  chosen  to  approximate 
typical  results  from  printing  with  230  or  325  mesh  screens: 
a  single  printed  layer  is  usually  15-20  |xm,  while  a  double 
printed  layer  (dried  between  printings)  is  30 — 40  |xm  thick. 

2.1.  Resistance  loss  calculation 


The  overall  resistance  of  an  array  of  identical  cells,  illus¬ 
trated  in  Fig.  1,  was  obtained  by  calculating  the  resistance  of 
an  individual  repeat  period  ( A)  and  multiplying  by  the  num¬ 
ber  of  cells.  Component  resistances  R  were  calculated,  for  the 
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Fig.  1.  Cross  sectional  schematics  of  the  SS-SOFC.  Part  (a)  shows  an  indi¬ 
vidual  cell  with  current  flows  indicated.  Part  (b)  shows  an  overall  view  with 
dimensions  labeled.  Part  (c)  shows  a  closeup  view  of  the  shunting  current 
calculation.  Note:  drawings  are  not  to  scale. 


most  part,  based  on  R  =  p/patt/4,  where  p  is  the  component 
resistivity,  /path  is  the  length  of  the  current  path,  and  A  is  the 
cross  sectional  area.  The  exception  to  this  was  the  cell  resis¬ 
tance  /('cell,  in  the  portion  of  the  device  in  Fig.  la,  where 
the  calculation  is  complicated  by  the  variation  in  cell  current 
density  with  lateral  position  (x).  The  total  resistance  for  each 
repeat  unit  (of  length  A)  also  included  electrode  length  in  the 
gap  between  the  active  cell  area  and  the  interconnect,  such 
that: 


Ra  =  //cell  +  - 
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where  the  dimensions  are  as  given  in  Fig.  lb,  pca  is  the  cath¬ 
ode  resistivity,  pan  is  the  anode  resistivity,  pi  is  the  inter¬ 
connect  resistivity,  and  w  is  the  cell  width  (into  the  plane  of 
Fig.  1).  The  interconnect  length  was  assumed  small  and  the 
electrode  resistivities  were  assumed  small  enough  that  the 
current  density  across  the  interconnect  was  uniform,  allow¬ 
ing  the  use  of  the  simple  interconnect  term  in  Eq.  (1).  The 
anode  and  cathode  gap  resistance  terms  were  approximated 
simply  by  including  half  of  the  interconnect  length  in  the 
electrode  current  paths.  This  does  not  introduce  substantial 
error,  given  that  this  resistance  is  small  relative  to  the  other 
terms  in  Eq.  (1).  Maximum  power  density  (MPD)  was  cal¬ 
culated  using  (MPD)  =  0.552 / RaAw  where  an  open  circuit 
voltage  of  1.1  V  was  assumed. 

The  calculation  of  //cell  accounted  for  the  anode  cur¬ 
rent  decrease  and  the  cathode  current  increase  due  to  current 
traversing  the  electrolyte  (Fig.  la).  The  equations  describing 
the  current  flow  and  voltage  drops  are  given  in  Eqs.  (2)— (4). 
The  parameters  are  defined  in  Appendix  A. 


d/an  —  d/ca  —  Jc\  /.(  i  i  |  dx 

Voc  -  (Vca  -  Van) 

= - Tcell  d-v 

RcM  +  //an/el  +  Rcl 


(2) 


d Van  =  -4n R  =  ~/an  ,  —  dx  (3) 

TcF.II.  fan 

dVca  =  /ca  R  =  -/cay— ^ — - dx  (4) 

TcF.l  t.  fca 

Integration  of  these  equations  was  performed  in  Matlab  [5] 
by  a  simple  iterative  technique.  The  initial  calculation  began 
at  x  =  0  with  the  values  /ca(x  =  0)  and  Van(x  =  0)  set  to  0,  cell 
current  (/cell)  set  to  /an  (x  =  0),  and  a  trial  value  of  Vca(x  =  0) 
chosen.  Eq.  (2)  gives  d/an  =  — d/ca,  calculated  as  the  current 
flow  across  the  electrolyte  between  x  and  x  +  dx.  For  simplic¬ 
ity,  the  polarization  resistances  /ZCa/ei  and  //an/e l  were  assumed 
to  be  ohmic.  Eqs.  (3)  and  (4)  give  the  voltage  drops  from  x 
to  x  +  dx  for  the  anode  dVan  and  cathode  dVca,  respectively. 
Vca(x=0)  was  iteratively  varied  using  a  bisection  algorithm 
until  the  boundary  condition  that  /an(x  =  Tcell)  =  0  was  sat¬ 
isfied.  Calculations  with  different  length  increments  dx  were 
done  and  the  results  compared  to  ensure  convergence  to  the 
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same  value.  The  resulting  values  were  used  to  give  the  re¬ 
sistance  flcFT  i  =  V™(x  =  LcFn)/IcPU.  in  Eq.  (1).  The  cell’s 
operating  voltage  is  defined  as  the  difference  between  the 
cathode  potential  at  x  =  Lcell  and  the  anode  potential  at  x  =  0, 
Van(x  =  0)  =  0,  such  that  Vtf.i.t.  =  VCd(x  =  Lcell)- 

As  a  check  of  the  present  numerical  calculation,  total  re¬ 
sistance  values  were  compared  with  those  obtained  from  the 
analytical  expression  given  in  Ref.  [3],  and  the  agreement 
was  very  good. 


2.2.  Shunting  current  calculation 


This  section  describes  the  calculation  of  the  current  flow¬ 
ing  between  cells  through  a  mildly  conductive  support,  as  il¬ 
lustrated  in  Fig.  lc.  The  goal  was  to  provide  a  more  accurate 
result  than  in  previous  work  [4],  where  the  crude  assumption 
of  a  uniform  field  across  the  support  was  employed.  Fig.  lc 
shows  a  cross  sectional  schematic  of  two  cells  along  with 
the  assumed  potential  versus  lateral  position  at  the  surface 
of  the  support  (y  =  0).  Since  a  typical  Ni-YSZ  anode  is  quite 
conductive,  the  potential  gradient  is  entirely  across  the  gap 
between  anodes: 

V(x)  =  -—  for  -L<x<-1^^  (5) 

2  ~  ~  2 

T //  x  r  ^G,an  ^  ,  ^G,an 

V(x)  =  - x  for - <  x  <  -  (6) 

^G,an  2  2 

V(x)  =  V°  for  lf^<x<L  (7) 

where  V0  is  the  cell  operating  voltage.  Since  the  potential  was 
assumed  not  to  vary  in  the  z  direction  (perpendicular  to  the 
plane  of  the  cross  section  pictured  in  Fig.  lc),  it  was  approx¬ 
imated  by  a  two-dimensional  solution  to  Faplace’s  equation 
for  a  rectangular  area  [6]: 

A0  ^ — -\ 

V(x,  y)  =s  —  +  y  (An  cos  finx  +  Bn  sin  finx) 

~  n 

X(  Cnep"y  +  Dne~^y)  (8) 


where  the  constants  Aa,An,Bn,  and  fin  can  be  calculated  using 
the  case  at  y  =  0  when  Yn  =  Cne^n}'  +  Z)„e-^"-v  =  1: 
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In  order  to  solve  for  the  constants  Cn  and  D„,  two  boundary 


conditions  were  used: 

II 

o' 

II 

(13) 

at  y  =  -ts,  ——  =  0 
dy 

(14) 

The  first  boundary  condition  assures  that  potential  profile 
will  match  Eqs.  (5) — (7)  at  y  =  ( ),  while  the  latter  condition 
expresses  that  no  current  can  flow  in  the  y  direction  at  the  free 
surface  of  the  support.  The  calculation  was  carried  out  over 
four  complete  cells,  and  the  potential  variation  taken  between 
the  central  two  cells  in  order  to  eliminate  edge  effects. 

The  expression  for  the  potential  profile  was  then  differen¬ 
tiated,  according  to  the  expression 


Ps 

in  order  to  find  the  approximate  current  density  between 
the  two  cells  [7],  The  expression  is  applicable  to  weakly- 
conductive  ohmic  materials  with  isotropic  dielectric  prop¬ 
erties,  where  localized  charges  do  not  affect  the  potential 
gradient. 


3.  Results 

3.1.  Resistance  loss  calculation 

In  this  section,  the  calculated  voltage  and  current  distri¬ 
butions  across  SS-SOFCs  are  described  for  two  categories 
of  cells  where:  (a)  the  cathode  and  anode  sheet  resistances 
are  approximately  equal,  and  (b)  one  electrode  has  a  sig¬ 
nificantly  smaller  sheet  resistance.  More  attention  is  given 
to  this  latter  case,  since  it  is  a  more  common  situation  in 
actual  SOFCs,  in  which  the  Ni-based  anode  is  more  con¬ 
ductive  than  the  perovskite  cathode.  For  Ni-based  anodes  a 
sheet  resistance  of  0.5  S2/D  is  used  based  on  a  thickness  of 
20  |xm  and  a  resistivity  of  ~  1  x  10~3  Q  cm  [8],  as  previously 
reported  for  porous  Ni-(Zr02)o.92(Y203)o.o8  (Ni-YSZ)  in 
a  55:45  volume  ratio  [4],  For  perovskite  cathode  materials 
such  as  (La,Sr)Mn03  (FSM),  a  sheet  resistance  of  10  Q/O  is 
used,  based  on  a  resistivity  of  «0.02  £2  cm  and  a  20  p,m  layer 
thickness.  For  simplicity,  the  individual  cells  are  assumed  to 
follow  ohmic  behavior,  and  the  cell  area-specific  resistance 
Ras  =  Rca/ei  +  Ran/el  +  Re l  value  was  0.5  £2  cm2  unless  other¬ 
wise  specified. 

3.1.1.  Symmetrical  electrode  cell  current  variations 

This  section  describes  calculations  done  for  the  case  where 
the  electrode  sheet  resistances  are  equal.  All  electrode  po¬ 
tential  profiles  and  electrolyte  current  densities  were  cal¬ 
culated  for  cells  operating  at  maximum  power,  i.e.  with 
Vcei.t.  =  VcAx= LCfu.)  =  0-55  V.  Fig.  2  shows  the  potential 
(a)  and  electrolyte  current  density  (b)  versus  position  for 
2  mm  long  cells,  for  sheet  resistances  of  0.5  or  10  S2/D.  For 
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the  lower  sheet  resistance,  the  potential  variation  across  each 
electrode  was  very  small,  with  a  nearly  constant  electrolyte 
current  density,  i.e.  there  were  effectively  no  electrode  re¬ 
sistance  losses.  For  the  higher  sheet  resistance,  however,  the 
potential  drop  was  about  0.15  V  across  the  2  mm  long  cell 
for  each  electrode.  The  potential  gradient  is  high  in  the  an¬ 
ode  and  low  in  the  cathode  near  x  =  0,  due  to  the  much  higher 
current  in  the  anode.  The  situation  is  reversed  near  x  =  Tcell, 
where  most  of  the  current  is  in  the  cathode.  The  electrolyte 
current  density  dropped  to  a  minimum  at  the  midpoint  of 
the  cell,  where  the  electrode  potential  difference  deviates  the 
least  from  OCV. 

Fig.  3  shows  the  results  of  a  similar  calculation  done  with 
10  mm  long  cells.  The  results  are  qualitatively  similar,  but  the 
larger  electrode  lengths  lead  to  larger  potential  variations. 
The  potential  variation  was  especially  large  for  the  higher 
sheet  resistance,  «^0.4  V  across  each  electrode,  correspond¬ 
ing  to  most  of  the  cell  potential  drop  of  0.55  V.  The  elec¬ 
trolyte  current  densities  displayed  more  variation  and  were 
lower  than  in  Fig.  2,  especially  for  the  larger  sheet  resistance 
where  the  current  density  dropped  to  near  zero  at  the  center 
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Fig.  2.  Anode  and  cathode  potential  (a)  and  electrolyte  current  density  (b) 
vs.  position  x  for  2-mm-long  symmetric-electrode  cells,  for  sheet  resistances 
of  0.5  or  lOfi/IH. 
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Fig.  3.  Anode  and  cathode  potential  (a)  and  electrolyte  current  density  (b)  vs. 
position  x  for  10-mm-long  symmetric-electrode  cells,  for  sheet  resistances 
of  0.5  or  10Q/D. 


of  the  cell.  These  results  show  that  higher  sheet  resistances 
require  the  use  of  smaller  cell  lengths  to  minimize  ohmic 
losses. 


3.1.2.  Asymmetrical  electrode  cell  current  variations 

These  calculations  were  done  for  sheet  resistance  values  of 
0.5  £2/D  for  the  anode  (e.g.  Ni-YSZ),  with  10  £!/□  or  30  £!/□ 
for  the  cathode,  the  latter  expected  for  »20  p,m  thick  LSM 
orLao.6Sro.4Coo.2Feo.gO3  (LSCF)  layers  taking  porosity  into 
consideration  [9].  Fig.  4  illustrates  the  electrode  potential  (a) 
and  electrolyte  current  density  (b)  versus  position  across  a  cell 
at  maximum  power  output.  The  anode  potential  was  essen¬ 
tially  constant.  The  cathode  potential  did  vary  substantially 
with  position.  For  the  cathode  sheet  resistance  of  30  £!/□,  the 
cathode  potential  dropped  ~0.3  V,  as  compared  to  ~0. 15  V 
for  the  lower  cathode  sheet  resistance.  As  shown  in  Fig.  4b, 
the  electrolyte  current  for  the  high  sheet  resistance  case  was 
smaller  than  the  low  resistance  case  near  j-Omm,  as  ex¬ 
pected  because  the  voltage  was  closer  to  Vqc  (Fig.  4a).  The 
electrolyte  current  variation  was  larger  for  the  higher  cathode 
sheet  resistance. 
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Fig.  4.  Anode  and  cathode  potential  (a)  and  electrolyte  current  density  (b) 
vs.  position  x  for  2-mm-long  asymmetric-electrode  cells,  for  cathode  sheet 
resistances  of  10  or  30  Q/D. 


For  comparison,  the  same  calculation  was  carried  out  for 
Lcell  =  10  mm.  Fig.  5  shows  the  electrode  potentials  (a)  and 
electrolyte  current  density  (b)  versus  length.  The  overall  de¬ 
pendence  is  similar  to  that  in  Fig.  4,  although  the  potential 
gradients  and  the  increase  in  current  density  at  large  x  are 
more  pronounced.  Most  of  the  current  travels  through  the  an¬ 
ode  and  crosses  the  electrolyte  at  large  x  in  order  to  avoid  a 
long  current  path  across  the  more  resistive  cathode.  For  the 
lower  cathode  sheet  resistance,  the  electrolyte  current  initially 
decreases  with  increasing  x,  shows  a  weak  minimum  near 
x  =  3  mm,  and  then  increases.  The  initial  decrease  near  x  =  0 
occurs  because  of  a  slight  decrease  in  the  local  cell  voltage, 
i.e.  a  decrease  in  anode  potential  while  the  cathode  poten¬ 
tial  remains  constant  (Fig.  5a).  The  anode  potential  change 
occurs,  despite  the  relatively  low  sheet  resistance  of  the  an¬ 
ode,  because  almost  all  of  the  cell  current  is  in  the  anode  at 
x  =  0.  There  is  little  corresponding  cathode  potential  change 
because  little  current  is  in  the  cathode.  As  x  increases,  more 
of  the  current  transfers  to  the  cathode,  such  that  the  cathode 
potential  decreases  rapidly.  A  comparison  of  Figs.  4b  and  5b 
shows  that  the  high  cathode  sheet  resistance  prevents  substan¬ 
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Fig.  5.  Anode  and  cathode  potential  (a)  and  electrolyte  current  density  (b) 
vs.  position  x  for  10-mm-long  asymmetric-electrode  cells,  for  cathode  sheet 
resistances  of  10  or  30 


tial  current  flow  except  within  ~2  mm  of  the  interconnect  for 
both  the  2  and  10  mm  cells. 

Fig.  6  summarizes  results  on  electrolyte  current  density 
versus  position  for  total  cell  lengths  ranging  from  1  to  10  mm. 
It  is  clear  from  this  figure  that  smaller  cell  lengths  resulted  in 
more  uniform  current  distributions  and  higher  average  elec- 


Distance  Along  Cell  Length  (mm) 


Fig.  6.  Current  density  through  the  electrolyte  for  cell  lengths  from 
1-10  mm.  Calculations  were  based  on  cell  operation  at  maximum  power 
density,  cathode  sheet  resistance  of  10  Q/d,  and  Ras  =  0.5  £2/d. 
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Active  Cell  Length--LCELL  (mm) 


Fig.  7.  Maximum  power  densities  vs.  cell  length  for  various  area  specific 
resistances  and  cell  lengths.  Only  active  cell  area  was  used  and  interconnect 
losses  were  not  included  in  this  power  density  calculation.  The  cathode  sheet 
resistance  was  10  S2/D. 

trolyte  current  densities,  primarily  due  to  the  shorter  electrode 
lateral  current  paths. 


3.1.3.  Performance  of  asymmetrical  cells 

In  this  section,  we  illustrate  the  impact  of  electrode  resis¬ 
tance  losses  on  overall  performance.  Fig.  7  shows  the  max¬ 
imum  power  density  (MPD)  versus  active  cell  length  Lcell 
for  an  anode  sheet  resistance  of  0.5  S2/EU,  a  cathode  sheet  re¬ 
sistance  of  10  S2/D,  and  Ras  =  1.0,  0.5,  0.3,  and  0.2  £2  cm2. 
In  this  figure,  power  densities  were  calculated  using  only 
the  active  cell  area.  The  greatest  maximum  power  densities 
occurred  at  the  smallest  cell  lengths,  as  expected  given  the  re¬ 
sults  in  the  previous  sections.  Cell  lengths  of  10  mm  resulted 
in  much  smaller  power  densities. 

The  calculation  was  extended  to  include  losses  due  to  the 
interconnect  and  inactive  portions  of  the  electrodes  (Eq.  (1)), 
and  account  for  the  full  surface  area  of  the  array.  Table  1 
gives  the  values  used  in  the  calculations.  Referring  to  Fig.  1, 
the  surface  can  be  divided  into  active  cell  length  (Tcell) 
and  interconnect/gap  length  (L\  +  L(',an  +  L(,ca).  We  have  as¬ 
sumed  interconnect/gap  lengths  that  are  fixed  by  processing 
considerations.  For  example,  in  screen  printing,  there  are  lim¬ 
itations  on  resolution  and  print-to-print  alignment.  Two  cases 
are  considered  here.  First,  lengths  were  assumed  that  are  rel¬ 
atively  easy  to  achieve,  i.e.  Lo.an  =  Lc,.cd  =  L\  =  0.25  mm.  Sec¬ 
ond,  we  chose  lengths  near  the  limits  of  screen  printing,  i.e. 
Eg, an  =  Tc.ca  =  U  =  0.1  mm. 

Fig.  8  shows  MPD  curves  for  the  two  cases,  each  showing 
an  optimal  MPD  for  each  Ras  value.  Indicated  in  the  figure 


Table  1 


Parameters  used  in  maximum  power  density  calculations 


Parameter 

Value 

Anode  sheet  resistance 

0.5  £2/n 

Cathode  sheet  resistance 

lOQ/n 

Interconnect  resistivity 

1  Q  cm 

Interconnect  thickness 

20  |xm 

are  the  values  of  Lcell  corresponding  to  the  optimal  MPDs. 
Also  shown  on  these  figures  is  the  fraction  of  the  surface  that 
is  active  cell  area.  The  MPD  occurs  at  active  area  fractions 
from  ~0.5  to  0.8.  The  maximum  is  a  result  of  the  electrode  re¬ 
sistance  effect,  shown  in  Fig.  7,  combined  with  the  decreasing 
fraction  of  active  cell  area  with  decreasing  Lcell-  Maximum 
power  densities  were  substantially  lower  than  the  maximum 
achievable  shown  in  Fig.  7.  Decreasing  the  gap/interconnect 
lengths  allowed  narrower  cells  while  maintaining  a  reason¬ 
ably  high  active  cell  fractional  area.  Thus,  the  maximum 
power  densities  were  higher  and  optimal  cell  lengths  smaller 
for  the  smaller  gap/interconnect  lengths  shown  in  Fig.  8b. 
Note  that  for  the  smaller  gap/interconnect  lengths,  the  elec¬ 
trode  lengths  are  decreased,  decreasing  that  resistance,  but 
the  narrower  interconnect  has  a  higher  resistance. 

Besides  the  cell  length  and  Ras,  the  other  key  factor 
affecting  cell  performance  is  the  cathode  sheet  resistance. 
Since  typical  cathode  material  resistivities  fall  within  a  nar¬ 
row  range,  the  sheet  resistance  will  be  primarily  varied  via 
cathode  thickness.  Fig.  9  shows  the  MPD  versus  cell  length 
for  different  cathode  sheet  resistances.  Decreasing  the  sheet 
resistance  allows  good  power  densities  to  be  maintained  at 
higher  cell  lengths.  For  example,  for  a  10  mm  cell  with  a 
10  Q/D  cathode  sheet  resistance  (~20  p,m  thick),  the  MPD  is 


(b)  Active  Cell  Length~LCELL  (mm) 


Fig.  8.  Maximum  power  density,  calculated  including  cell  and  interconnect 
area,  for  (a)  Lo.an  =  i-G.ca  =  Li  =  0.25  mm  and  (b)  Lo.an  =  f-G.ca  =  U  =0.1  mm 
using  different  values  of  R\s  ■  Also  shown  is  the  percent  active  area  for  each 
cell  length. 


T.S.  Lai,  S.A.  Barnett  /  Journal  of  Power  Sources  147  (2005)  85-94 


91 


Fig.  9.  Maximum  power  density,  calculated  including  cell  and  interconnect 
area,  as  a  function  of  active  cell  length  for  various  cathode  sheet  resistances 
and  R\s  =  0-5  £2  cm2.  The  present  screen  printed  cathode  sheet  resistance  is 
«10  Q/D. 

only  118  mW  cm-2.  Decreasing  the  sheet  resistance  to  1  S2/D 
(~200  |jLm  thick),  the  MPD  is  increased  to  307  mW  cm-2.  It 
should  be  noted  that  MPDs  at  a  length  of  10  mm  are  less  than 
those  of  narrower  cells,  and  that  a  processing  technique  other 
than  screen  printing  would  be  required  to  make  such  thick 
cathode  layers. 

Resistance  losses  in  the  interconnect  region  can  also  in¬ 
fluence  overall  cell  performance.  Because  of  the  relatively 
high  resistivity  of  typical  ceramic  interconnect  materials,  this 
can  be  an  important  effect.  To  demonstrate  the  effect  of  the 
interconnect  length  on  array  performance,  the  array’s  opti¬ 
mum  MPD  and  the  optimum  active  cell  length  are  plotted 
versus  interconnect  length  in  Fig.  10.  Ras  was  assumed  to 
be  0.5  £2  cm2  and  the  electrode  gap  lengths  were  kept  con¬ 
stant  at  Lq  ca  =  /-G  an  =  0.25  mm.  With  those  parameters,  an 
interconnect  length  of  0. 1  mm  and  an  associated  active  cell 
length  of  1.35  mm  gave  the  highest  MPD.  Larger  intercon¬ 
nect  lengths  resulted  in  lower  MPDs  due  to  the  smaller  frac¬ 
tion  of  active  cell  area.  Decreasing  the  length  below  0. 1  mm 


Fig.  10.  Maximum  power  density  and  optimum  active  cell  length  vs.  inter¬ 
connect  length  for  interconnect  thickness  of  20  p.m  and  resistivity  of  1  Q  cm. 


Fig.  1 1 .  Equipotential  plots  for  gap  lengths  (a)  0. 1  mm  and  (b)  1  mm.  Gaps 
are  bounded  by  arrows. 


also  lowered  MPDs  because  of  the  increasing  interconnect 
resistance. 

3.2.  Shunting  current  calculation 

Fig.  11  shows  examples  of  the  calculated  equipotential 
lines  within  the  support  for  cells  with  Lcell  =  2  mm,  sup¬ 
port  thickness  fs  =  1  mm,  and  gap  lengths  Lc.an  =  0.1  (a)  and 
1  mm  (b).  From  Eq.  (11),  current  flow  is  perpendicular  to 
the  equipotential  lines,  with  magnitude  proportional  to  the 
line  density.  Directly  below  each  electrode,  there  were  few 
equipotential  lines  and  hence  little  current  flow,  as  expected 
since  the  electrodes  are  at  constant  potential.  Well  below  the 
electrode  gap,  the  lines  were  generally  uniformly  spaced,  in¬ 
dicating  uniform  current  flow.  Near  the  electrode  gaps,  the 
lines  were  more  tightly  spaced,  indicating  a  higher  current 
density.  For  the  smaller  gap,  the  lines  were  more  tightly 
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bunched,  indicating  a  higher  field  and  current  density  local¬ 
ized  in  a  smaller  volume. 

In  order  to  provide  an  example  of  the  shunting  cur¬ 
rent  expected  in  a  typical  SS-SOFC,  we  have  assumed 
PSZ  supports.  Resistivity  values  for  dense  PSZ  range  from 
^1300-220  £2  cm  between  575  and  1000  °C,  respectively 
[10].  These  values  were  corrected  to  account  for  the  typical 
porosity  of  30-40  vol.%  [4].  Since  PSZ  is  ionically  conduc¬ 
tive  and  the  Ni-YSZ  electrodes  are  electronically  conductive, 
current  flow  through  the  support  requires  electrochemical  re¬ 
actions  at  each  of  the  interfaces  (Fig.  lc),  resulting  in  an  as¬ 
sociated  polarization  resistance.  This  resistance  is  neglected 
in  the  present  calculation  since  its  value  may  vary  greatly  de¬ 
pending  on  exact  electrode  composition  and  structure.  Thus, 
the  shunting  current  values  given  below  must  be  regarded  as 
upper  limits. 

The  net  substrate  current  /shunt was  obtained  by  integrating 
the  current  in  Eq.  (15)  over  the  thickness  of  the  support  at  the 
midpoint  between  the  two  cells  (x  =  0  in  Fig.  11): 


f  shunt  =  J(x  =  0jd-dy  = 


z,y 


z,y 


-V  V(x  =  0) 
PS 

r  dv 

dx 


djdy 

dy  (16) 


x=0 


Note  that  since  d  V7dy  =  d Vldz  =  0  at  the  midpoint  by  symme¬ 
try,  VV=dV7dv.  We  have  plotted  /shuntPs/w  in  Fig-  12  as  a 
function  of  Lq  w  for  three  different  values  of  Lcell  :  F  2,  and 
5  mm  and  a  support  thickness  of  1  mm.  The  shunting  current 
increased  gradually  with  decreasing  gap  length.  This  can  be 
explained  from  the  equipotential  plots  in  Fig.  11.  While  the 
equipotential  lines  far  from  the  gap  changed  little  with  gap 
length,  the  lines  became  very  closely  spaced  within  the  nar¬ 
rower  gap,  indicating  a  very  high  local  current  density  that 
increased  the  overall  current.  The  cell  length  had  little  effect 
on  the  shunting  current.  This  can  be  understood  by  noting 
that  the  equipotential  lines  in  the  gap  region  varied  little  with 
cell  length,  and  this  region  was  the  primary  factor  limiting 
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Fig.  12.  Shunting  current  vs.  gap  length  for  various  active  cell  lengths. 


Active  Cell  Length~LCELL  (mm) 


Fig.  13.  Ratio  of  shunting  current  to  cell  current  plotted  vs.  active  cell  length 
for  gap  lengths  of  0. 1  and  0.25  mm.  Cell  currents  were  obtained  from  the 
data  plotted  in  Fig.  7,  assuming  7?as  =0.5  £2  cm2. 


the  shunting  current  (in  the  active  cell  length  the  conductive 
Ni-YSZ  anode  provided  an  easy  short-circuit  current  path). 
The  present  results  are  for  one  specific  support  thickness;  in 
general,  the  shunting  current  will  increase  with  increasing 
support  thickness. 

Fig.  13  shows  the  ratio  of  7shUnt  to  /cell  as  a  function 
of  Lcell,  assuming  an  operating  temperature  of  800  °C. 
Le  an  was  taken  to  be  0.25  mm  and  cell  performance  data 
were  the  same  as  the  R,\s  =  0.5  !T2  cm2  series  in  Fig.  7.  At 
/-cell  ~6  mm,  the  ratio  goes  through  a  minimum  value 
of  about  0.026  due  to  the  increasing  cell  current.  Be¬ 
low  Lcell  mm,  /shunt//cELL  >0.05.  As  Lckll  approaches 
10  mm,  the  ratio  rises  slightly  because  the  long  current  path 
has  a  detrimental  effect  on  cell  performance. 

Fig.  14  shows  the  predicted  shunting  current  assuming 
LG,an  =  0.4mm,  Lcell  =  1.8  mm,  and  w  =  8.5mm.  The  cur¬ 
rent  increases  with  increasing  cell  operating  temperature. 
The  highest  value  was  8.2  mA  at  1000  °C,  whereas  the  low- 
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Fig.  14.  Shunting  current  vs.  temperature  comparing  the  present  calcula¬ 
tion,  a  previous  approximate  calculation  assuming  a  uniform  field  across 
the  support  [4],  and  an  experimental  measurement.  For  each  of  these, 
Lc,an  =  0.4  mm,  Lq ell  =  1.8  mm,  and  w  =  8.5  mm. 
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est  value  was  1.4  mA  at  600  °C.  At  the  expected  operating 
temperature  of  800  °C  (a  resistivity  of  752  cm  after  correc¬ 
tion  for  porosity),  the  shunting  current  was  calculated  to  be 
4.3  mA,  translating  to  at  most  a  5-6%  degradation  in  MPD 
when  considering  the  performance  data  given  in  Fig.  7. 

An  experimental  check  of  the  calculation  results  was  ob¬ 
tained  by  making  resistance  measurements  (in  humidified  hy¬ 
drogen)  across  an  array  of  four  Ni-YSZ  anodes  printed  onto 
the  surface  of  a  porous  PSZ  substrate.  The  dimensions  were 
the  same  as  those  used  in  the  calculation  in  Fig.  14.  The  mea¬ 
sured  current,  shown  in  Fig.  14,  increased  with  increasing 
temperature,  like  the  calculation  result,  but  at  a  faster  rate. 
While  the  reason  for  the  difference  in  slope  is  not  known, 
there  was  reasonable  agreement  between  the  magnitudes  of 
the  measured  and  predicted  shunting  currents.  Note  that  the 
measured  values  are  lower  than  the  predicted  values  at  lower 
temperatures;  this  may  be  due  to  the  effect  of  polarization 
resistance,  which  was  not  considered  in  the  calculation.  Also 
plotted  are  shunting  current  values  calculated  assuming  a  uni¬ 
form  electric  field  across  the  array  repeat  period  A  and  an  op¬ 
erating  voltage  of  0.55  V,  as  mentioned  in  our  previous  paper 
[4].  This  method  produced  shunting  current  values  at  least 
a  factor  of  10  less  than  those  calculated  using  the  present 
method.  These  results  illustrate  the  importance  of  operat¬ 
ing  temperature  on  the  shunting  current,  originating  from  the 
temperature  dependence  of  the  PSZ  resistivity.  Lower  tem¬ 
peratures  are  clearly  desirable  to  minimize  shunting  current, 
assuming  that  the  cell  current  does  not  decrease  too  much  at 
the  lower  temperature. 

4.  Discussion 

The  above  results  show  how  the  geometry  of  segmented- 
in-series  SOFCs  determine  their  performance.  For  cells  with 
low  area  specific  resistance,  relatively  small  cell  lengths  and 
cathode  sheet  resistances  are  required  to  achieve  high  power 
densities,  in  agreement  with  prior  predictions  [4].  Otherwise, 
the  electrolyte  current  is  non-uniform  across  the  length  of  the 
cells,  with  much  of  the  area  not  contributing  significantly  to 
the  cell  current.  The  above  results  predict  an  optimum  cell 
length  by  accounting  for  the  interconnect  area.  Assuming 
typical  SOFC  materials  properties  and  corresponding  area- 
specific  resistance  values  of  «0.3-0.5  cm2,  along  with  ge¬ 
ometries  and  layer  thickness  characteristic  of  screen  print¬ 
ing,  optimal  power  densities  were  achieved  for  cell  lengths 
of  1-2  mm.  For  the  10-mm-long  cells  that  have  traditionally 
been  used  in  segmented-in-series  SOFCs,  power  densities  are 
much  lower  unless  very  thick  (>200  p,m)  cathodes  or  more 
conductive  (^>100S  cm-1)  cathode  materials  are  employed. 

Interconnect  resistances  for  segmented-in-series  SOFCs 
did  not  become  a  limiting  factor,  despite  the  relatively  low 
conductivity  («H  S  cm-1)  assumed  for  typical  ceramic  inter¬ 
connects,  unless  the  interconnect  length  was  <0.1  mm  (<5% 
of  the  cell  area).  This  is  primarily  due  to  the  segmented-in- 
series  geometry,  in  which  the  thin  («20  |Jim)  interconnect  is 


sandwiched  between  more  conductive  anode  and  cathode  lay¬ 
ers.  Since  standard  materials  such  as  LaCr03  sinter  poorly, 
thin  and  dense  ceramic  interconnect  layers  are  difficult  to  pro¬ 
duce  [1 1].  In  addition,  the  interconnect  should  be  co-sintered 
with  the  support,  anode,  and  electrolyte  layer,  providing  fur¬ 
ther  constraints  on  sintering  conditions  to  avoid  possible  reac¬ 
tions  that  could  form  resistive  interfacial  phases.  Nonetheless, 
dense  LaCrC>3  -based  interconnect  layers  have  been  reported 
using  sintering  aids  [12]. 

The  effects  of  shunting  current  were  evaluated  for  the 
case  of  zirconia-based  supports,  which  are  convenient  for 
SS-SOFCs  because  of  the  good  chemical  stability  and  ther¬ 
mal  expansion  match  with  SOFC  materials.  For  reasonable 
area-specific  resistance,  cell  length,  and  operating  tempera¬ 
ture  <800  °C,  the  shunting  current  will  not  substantially  af¬ 
fect  overall  power  density.  However,  shunting  current  may 
become  important  for  higher  temperatures  and  for  cell  re¬ 
peat  periods  <1.2  mm,  where  the  shunting  current  com¬ 
prises  >  5%  of  the  cell  current. 

While  it  is  inevitable  that  comparisons  will  be  made  be¬ 
tween  area-specific  power  densities  for  different  stack  de¬ 
signs,  volumetric  power  density  and  cost  per  kW  are  the 
more  relevant  quantities.  For  this,  the  different  stack  geome¬ 
tries  must  be  considered.  While  these  considerations  are  be¬ 
yond  the  scope  of  the  present  work,  a  few  qualitative  re¬ 
marks  can  be  made  regarding  how  to  meaningfully  compare 
the  present  power  densities  with  planar  SOFCs.  First,  SS- 
SOFC  power  densities  calculated  here  include  interconnect 
and  current  collection  losses,  such  that  they  should  be  com¬ 
pared  with  planar  SOFC  stuck  results,  rather  than  single  cells. 
Note  that  the  power  densities  predicted  above  for  SS-SOFCs 
are  generally  lower  than  for  planar  single  cells,  which  are 
typically  >  1  W  cm-2  at  800  °C,  but  are  similar  to  those  for 
planar  SOFC  stacks  (typically  rs  0.5  W  cm-2).  Second,  the 
power  densities  in  Fig.  8  include  both  cell  and  interconnect 
areas  in  the  calculation.  In  contrast,  power  density  calcula¬ 
tions  for  planar  stacks  generally  do  not  include  interconnect 
area.  Thus,  a  more  appropriate  comparison  would  be  to  use 
the  power  density  calculated  using  only  active  cell  area,  e.g. 
Fig.  7,  evaluated  at  an  appropriate  cell  length. 

5.  Summary  and  conclusions 

Calculations  were  done  to  predict  (a)  the  electrical  be¬ 
havior  of  segmented-in-series  SOFCs  based  on  their  geome¬ 
try  and  components’  sheet  resistances,  and  (b)  the  shunting 
current  flowing  between  cells  through  the  porous  substrate. 
Calculations  were  made  using  array  geometries  reflecting  the 
capabilities  of  screen  printing  and  conductivities  typical  of 
SOFC  materials.  For  a  gap  length  of  0.25  mm  and  an  as¬ 
sumed  operating  temperature  of  800  °C,  the  shunting  current 
through  the  substrate  was  estimated  to  be  <6%  of  typical 
cell  currents.  For  area  specific  cell  resistances  of  0.2,  0.3, 
0.5,  and  1.0  cm2,  the  optimum  maximum  power  densities 
occurred  from  1 .05-1 .9  cm  for  Lo.an  =  Lo.ca  =  k\  =  0.25  mm 
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and  0.8-1 .45  cm  for  Lq  w  =  Lq  ca  =  L\  =  0. 1  mm.  The  calcu¬ 
lations  can  be  used  to  design  and  optimize  SS-SOFCs  for 
a  given  materials  set,  operating  condition,  and  processing 
technique.  For  typical  SOFC  materials  and  layer  thicknesses 
typically  obtained  by  screen  printing,  reducing  the  active  cell 
length  from  10  mm  to  1-3  mm  should  result  in  a  significant 
improvement  in  power  density. 
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Appendix  A.  Definitions  for  symbols 


Symbol 

Meaning 

d/an 

Change  in  anode  current  in  length  element  dr 

d/ca 

Change  in  cathode  current  in  length  element  dr 

dVan 

Voltage  drop  across  anode  in  length  element  dr 

dVca 

Voltage  drop  across  cathode  in  length  element  dr 

/an 

Current  flowing  through  anode 

/ca 

Current  flowing  through  cathode 

/cell 

Cell  current 

/shunt 

Shunting  current 

J el 

Flux  of  electrons  thrown  electrolyte 

A 

Repeat  period  of  array 

L 

Anode  length  plus  half  of  anode  gap  length 

/'CELL 

Cell  length 

Z<J,an 

Gap  length  between  anodes 

Z'G.ca 

Gap  length  between  cathodes 

/'I 

Length  of  interconnect 

Pan 

Resistivity  of  anode 

Pea 

Resistivity  of  cathode 

^an/el 

Resistance  of  anode/electrolyte  interlace 

Symbol 

Meaning 

Ras 

Area  specific  resistance  of  cell  (excluding  elec¬ 
trodes),  Rel  +  Ran/el  +  Rca/el 

Rcajel 

Resistance  of  cathode/electrolyte  interlace 

Re  1 

Resistance  of  electrolyte 

tan 

Thickness  of  anode 

tea 

Thickness  of  cathode 

tel 

Electrolyte  thickness 

1\ 

Interconnect  thickness 

ts 

Substrate  thickness 

Van 

Voltage  of  anode  (relative  to  0) 

Vca 

Voltage  of  cathode  (relative  to  0) 

Vo 

Cell  operating  voltage 

Foe 

Open  circuit  voltage  of  cell 

w 

Cell  width 

X 

Position  along  cell  length 

y 

Distance  from  substrate  surface 

z 

Direction  perpendicular  to  xy  plane 
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